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ABSTRACT 


This  effort  focuses  on  collecting  and  analyzing  regional  seismic  data  for  earthquakes  and  underground  explosions  to 
improve  seismic  event  characterization  capabilities  with  regard  to  monitoring  the  Comprehensive  Nuclear-Test-Ban 
Treaty  (CTBT).  Digital  data  have  now  been  acquired  and  processed  for  267  regional  recordings  of  160  underground 
nuclear  (UNE)  and  chemical  (UCE)  explosions.  The  explosion  data  set  includes  99  Nevada  Test  Site  (NTS)  UNE’s 
recorded  by  MNV  and  KNB;  6  Lop  Nor  UNE’s  recorded  by  NIL,  ZAL,  ULN,  WMQ,  BRVK  and  MAKZ;  22 
Semipalatinsk  (STS)  UNE’s  recorded  by  WMQ  and  BRVK;  6  Novaya  Zemlya  (NZ)  UNE’s  recorded  by  ARCES 
and  KEV;  and  one  UNE  at  the  Indian  Test  Site  recorded  by  NIL.  In  addition,  the  set  includes  20  Former  Soviet 
Union  (FSU)  peaceful  nuclear  explosions  (PNE)  recorded  by  BRVK  and  KEV;  and  6  STS  UCE’s  recorded  by  ZAL 
and  BRVK.  (Note:  not  all  of  these  explosions  have  recordings  by  all  stations  listed.)  Over  5300  regional  recordings, 
between  3  and  17  degrees,  by  existing  seismic  stations  of  the  International  Monitoring  System  (IMS)  are  also  used 
in  this  study.  The  regional  recordings  correspond  to  4173  presumed  earthquakes  above  mb  3.5  in  the  Reviewed 
Event  Bulletin  (REB)  that  were  processed  at  the  prototype  International  Data  Center  (pIDC).  Maximum  amplitude 
measurements  of  Pn,  Pg,  Sn  and  Lg  have  been  computed  in  several  frequency  bands  ranging  from  2  Hz  up  to  14  Hz 
(depending  on  the  Nyquist  frequency  of  a  given  station). 

The  data  for  the  presumed  earthquakes  have  been  used  to  estimate  attenuation  corrections  for  Pn/Sn  and  Pn/Lg 
amplitude  ratios  in  the  2^4,  4—6,  6-8,  and  8-10  Hz  passbands  for  tectonic  and  stable  region  types.  The  regional 
phase  amplitude  ratios  are  further  corrected  for  path  variations  using  a  Bayesian  Kriging  algorithm.  A  correction 
surface  is  determined  for  each  type  of  amplitude  ratio  at  each  station  as  an  optimal  linear  combination  of  existing 
amplitude-ratio  data  at  the  station,  giving  greater  weight  to  data  nearer  to  the  location  of  the  event  to  be  corrected.  A 
corresponding  uncertainty  surface  is  also  estimated  in  terms  of  the  residual  variance  of  the  data  and  a  calibration 
variance.  For  well-calibrated  locations,  the  correction  converges  to  the  mean  of  nearby  data  and  the  uncertainty 
converges  to  the  residual  variance.  For  locations  far  from  calibration  data,  the  correction  surface  converges  to  the 
worldwide  average,  with  larger  uncertainty  due  to  both  residual  and  calibration  variances. 

Using  these  correction  and  uncertainty  surfaces,  corrected  values  of  Pn/Smax(6-8  Hz)  are  then  used  to  define  a 
hypothesis  test  that  fixes  the  significance  level  with  respect  to  misclassifying  explosions.  A  score  is  computed  for 
each  event  such  that  events  with  negative  scores  are  not  screened  out  at  the  specified  significance  level.  The 
significance  level  used  here  has  been  set  such  that  none  of  the  available  underground  explosions  are  screened  out. 
Events  with  scores  greater  than  zero  are  screened  out  (i.e.,  are  considered  to  be  inconsistent  with  the  hypothesis  that 
the  event  belongs  to  the  explosion  population).  Screening  rates  for  earthquakes  above  mb  3.5  are  estimated  and 
improvements  in  event  characterization  performance  due  to  Kriging  are  demonstrated. 
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OBJECTIVE 


This  effort  focuses  on  collecting  and  analyzing  regional  seismic  data  for  underground  explosions  and  earthquakes  to 
improve  seismic  event  characterization  capabilities  with  regard  to  monitoring  the  Comprehensive  Nuclear-Test-Ban 
Treaty  (CTBT).  Key  aspects  of  this  effort  include:  (i)  acquiring  regional  seismic  recordings  for  underground 
explosions  and  earthquakes;  (ii)  processing  the  waveforms  to  measure  Pn,  Pg,  Sn  and  Lg  amplitudes  in  various 
frequency  bands;  (iii)  evaluating  optimum  processing  techniques  for  regional  amplitude  ratios;  (iv)  developing  a 
regional  event-screening  capability,  which  includes  establishing  necessary  region-specific  training  sets  and  distance 
corrections;  and  (v)  evaluating  optimum  screening  criteria  for  various  regions  and  comparing  explosion  and 
earthquake  characteristics  in  an  effort  to  understand  how  regional  event  characterization  parameters  and  screening 
criteria  can  be  transported  to  regions  for  which  no  explosion  data  exist. 

This  work  is  motivated  by  the  fact  that,  although  depth  and  Ms:mb  are  essential  methods  in  the  characterization  of 
seismic  events,  the  experimental  depth  and  Ms:mb  screening  criteria  being  tested  at  the  prototype  International  Data 
Center  (pIDC)  currently  screen  out  less  than  50%  of  the  presumed  earthquakes  above  mb  3.5  in  the  Reviewed  Event 
Bulletin  (REB)  (e.g.,  Fisk  et  ah,  1999b).  However,  there  are  a  significant  number  of  events  in  the  REB  with  regional 
seismic  data  that  can  provide  additional  event-screening  performance.  Many  studies  have  indicated  that  ratios  of 
high-frequency,  regional  seismic  phase  amplitude  ratios  can  provide  adequate  separation  between  underground 
explosions  and  earthquakes  to  aid  in  monitoring  the  CTBT  (see  Fisk  et  ah,  1996,  and  references  therein). 

We  describe  data  that  have  been  compiled  and  processed  for  267  regional  recordings  of  160  underground 
explosions.  These  explosions  are  compared  to  over  5300  regional  waveforms,  recorded  by  existing  seismic  stations 
of  the  International  Monitoring  System  (IMS),  for  4173  presumed  earthquakes,  above  mb  3.5,  processed  at  the 
pIDC.  After  applying  distance  corrections  that  depend  on  region  type  (tectonic  or  stable)  to  Pn/Sn  and  Pn/Lg  ratios, 
a  Bayesian  Kriging  algorithm  is  applied  at  each  station  to  further  correct  for  path  variations.  The  Kriging  algorithm 
also  provides  a  corresponding  uncertainty  surface.  To  characterize  the  events,  we  use  Pn/Smax  in  the  6-8  Hz  band, 
where  Smax  =  max(Sn.Lg),  which  has  been  corrected  for  distance  and  for  path  variations  using  the  Kriged  surface. 
We  apply  a  hypothesis  test,  which  accounts  for  the  varying  uncertainty  associated  with  the  Kriged  surface,  to  these 
data  sets  to  assess  whether  an  event  is  consistent  with  the  explosion  population  at  a  fixed  significance  level  with 
respect  to  incorrectly  screening  out  an  explosion.  A  score  is  computed  such  that  events  with  positive  scores  are 
screened  out,  i.e.,  are  inconsistent  with  the  hypothesis  that  the  event  belongs  to  the  explosion  population.  Events 
with  scores  less  than  or  equal  to  zero  are  not  screened  out.  At  the  0.005  significance  level,  none  of  the  available 
explosions  are  screened  out. 


RESEARCH  ACCOMPLISHED 


Regional  Explosion  Data 

We  have  compiled  267  regional  seismic  recordings  for  160  underground  nuclear  (UNE),  chemical  (UCE)  and 
peaceful  nuclear  (PNE)  explosions.  Table  1  summarizes  our  current  database  of  regional  explosion  recordings, 
which  represent  explosions  in  diverse  geological  regions  at  various  test  sites.  These  explosions  cover  a  magnitude 
range  of  about  ML  2.4  to  mb  6.2,  and  a  distance  range  of  about  2  to  17  degrees.  The  regional  measurements  for  the 
NTS  explosions  recorded  by  MNV  and  KNB  were  provided  by  LLNL  (Patton  and  Walter,  1994)  and  the  PNE  data 
recorded  by  BRVK  were  provided  by  Jack  Murphy  (1998). 

The  processed  data  consist  of  regional  phase  amplitudes,  Pn,  Pg,  Sn  and  Lg  in  the  2^1,  4—6,  6-8,  and  8-10  Hz 
bands.  Not  all  events  have  complete  sets  of  amplitude  ratios.  Six  of  the  UNE’s  conducted  on  Novaya  Zemlya  (NZ) 
were  recorded  by  station  KEV  and  three  were  recorded  by  ARCES.  Six  UNE’s  occurred  at  the  Lop  Nor  test  site,  one 
recorded  by  ZAL,  four  by  NIL,  two  by  ULN,  one  by  WMQ,  four  by  BRVK  and  four  by  MAKZ.  Semipalatinsk 
(STS)  provides  22  UNE’s,  21  recorded  by  WMQ  and  six  by  BRVK.  There  is  one  UNE  at  the  Indian  Test  Site 
recorded  by  NIL.  There  are  99  explosions  from  the  Nevada  Test  Site  (NTS),  84  of  which  have  recordings  at  both 
MNV  and  KNB,  nine  at  MNV  only,  and  six  at  KNB  only.  In  addition,  the  set  includes  20  Former  Soviet  Union 
(FSU)  PNE’s,  16  recorded  by  BRVK  and  six  by  KEV;  and  6  STS  UCE’s,  five  recorded  by  ZAL  and  six  by  BRVK. 
A  few  events  in  the  table  do  not  have  either  Pn/Sn(6-8  Hz)  or  Pn/Lg(6-8  Hz)  available,  and  therefore  not  are  used  in 
the  screening  test  described  below. 
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Table  1:  Summary  of  Regional  Underground  Nuclear  Explosion  Data 


Station 

Distance 

Magnitude 

Test  Site 

Explosion  Dates 

KEV 

9.4 

5.52-5.90 

NZ  UNE’s 

821011,841025,870802, 
880507,  881204,  901024 

3.1-19.5 

4.70-5.50 

FSU  PNE’s 

5 (820904-850718) 

ARCES 

9.9 

5.60-5.90 

NZ  UNE’s 

880507,  881204,  901024 

ZAL 

12.2 

4.71 

Lop  Nor  UNE 

960729 

6.2,  5.5,  5.5 

6.0,  6.0 

3.95,  n/a,  n/a, 

3.80,  3.67 

STS  UCE’s 

970803,970831,970928, 
980822,  990925 

NIL 

14.4-14.7 

5.73,  5.54 

5.69,  4.71 

Lop  Nor  UNE’s 

950515,950817, 

960608,  960729 

6.7 

5.00 

India  UNE 

980511 

ULN 

14.3 

5.73,5.54 

Lop  Nor  UNE’s 

950515,950817 

WMQ 

8.6 

4.6-6. 1 

STS  UNE’s 

21  (870606-891004) 

2.2 

4.7 

Lop  Nor  UNE 

880929 

BRVK 

16.8-16.9 

6.00,5.73 

5.54,  5.69 

Lop  Nor  UNE’s 

941007,  950515, 

950817,  960608 

6.2 

5.70-7.20 

STS  UNE’s 

6 (88023-881217) 

5. 8-6. 2 

3.95,  n/a,  n/a, 

3.80,  n/a,  3.67 

STS  UCE’s 

970803,970831,970928, 
980822,  980917,  990925 

7.2-17.2 

4.5-5. 8 

FSU  PNE’s 

16  (730815-880906) 

MAKZ 

6. 9-7. 2 

6.00,  5.73, 

5,69, 4.71 

Lop  Nor  UNE’s 

941007,  950515, 

960608,  960729 

MNV 

1. 7-2.2 

2.4-5 .5 

NTS  UNE’s 

93 (1979-1992) 

KNB 

2. 5-2.8 

2.6-5. 5 

NTS  UNE’s 

90  (1979-1992) 

Regional  Earthquake  Data 

Over  5300  regional  waveforms  that  were  recorded  by  52  IMS  seismic  stations  and  processed  at  the  pIDC  are  used  in 
this  study.  These  recordings  correspond  to  4173  distinct  REB  events  above  mb  3.5  that  are  presumed  to  be 
earthquakes.  The  restriction  to  mb  >  3.5  is  used  to  reduce  the  possibility  of  mining  blasts  contaminating  the 
earthquake  data  sets.  We  also  limit  the  events  to  those  with  distances  beyond  3  degrees  to  avoid  problems  of 
distinguishing  Pn  and  Pg,  and  to  those  less  than  17  degrees,  although  regional  phase  amplitudes  are  computed  at 
distances  up  to  20  degrees  at  the  pIDC.  In  addition,  we  limit  this  study  to  data  from  IMS  stations  that  have  at  least 
20  or  more  recordings  so  that  useful  training  sets  may  be  formed.  Figure  1  shows  the  locations  of  the  4173  presumed 
earthquakes  and  the  160  explosions  used  in  this  study.  The  locations  of  the  52  IMS  stations  and  4  non-IMS  stations 
(KEV,  WMQ,  BRVK  and  KNB)  are  included. 

Maximum  amplitudes  of  Pn,  Pg,  Sn  and  Lg  have  been  measured  at  the  pIDC  in  the  2—4,  4-6,  6-8,  8-10  and  10-12 
Hz  bands  (depending  on  the  Nyquist  frequency  of  the  instrument).  For  this  study,  we  use  Pn/Sn  and  Pn/Lg  in  the  6- 
8  Hz  bands  that  satisfy  the  signal-to-noise  criteria  of  Pn  signal  divided  by  pre-Pn  noise  greater  than  2.0  and  S  signal 
divided  by  pre-S  noise  greater  than  1.2. 

In  addition  to  the  pIDC  data,  the  earthquake  data  set  also  includes  60  recordings  by  MNV  and  75  by  KNB  of  Pn/Lg 
in  the  4-6,  6-8,  and  8-10  Hz  bands  for  earthquakes  near  NTS,  processed  by  LLNL  in  the  same  way  as  the  99 
UNE’s.  These  earthquakes  are  necessary  to  properly  calibrate  stations  MNV  and  KNB. 
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Figure  1.  Locations  of  160  nuclear  explosions  (triangles)  and  4173  REB  earthquakes  with  mb  3.5  and  above. 

Locations  of  52  IMS  stations  and  4  stations  with  explosion  recordings  (KEV,  WMQ,  BRVK,  and 
KNB)  are  also  shown. 

Distance  Corrections 

The  unequal  rates  of  attenuation  with  distance  between  regional  P  and  S  phases  necessitate  a  correction  for  distance 
if  amplitude  ratios  recorded  at  different  distances  are  to  be  compared  (see  Sereno,  1990,  Bottone  et  al.,  1997,  Fisk  et 
al.,  1998,  Fisk  et  al.,  1999a  and  Fisk  et  al.,  2000).  P/S  distance  dependence  is  modeled  by  a  three -parameter  equation 
of  the  form 

log(P/S)  =  a  +  blog  A  +  cA ,  (1) 

where  a,  b,  and  c  are  parameters  depending  on  the  particular  amplitude  ratio  (Pn/Lg  or  Pn/Sn),  frequency  band  and 
region  type.  A  is  the  epicentral  distance.  Following  Jenkins  et  al.  (1998),  we  divide  the  world  into  tectonic  and  stable 
region  types  and  estimate  the  three  coefficients  for  each  region  type.  We  then  correct  all  amplitude  ratios  at  a  given 
station  using  equation  (1)  with  coefficients  estimated  using  all  data  located  in  the  appropriate  region  type.  At  this 
point  we  do  not  make  station-dependent  corrections;  station  dependence  will  result  from  Rriging.  Figure  2  plots  the 
distance  dependence  of  all  REB  earthquakes  in  both  the  tectonic  and  stable  regions  for  Pn/Lg(6-8  Hz)  and  Pn/Sn(6- 
8  Hz). 

Bayesian  Kriging 

Corrections  of  regional  amplitude  ratios  using  worldwide  distance  corrections  do  not  depend  on  station  and  do  not 
take  into  account  any  variations  in  propagation  amplitude  that  may  depend  on  path  differences.  One  way  of  treating 
path  variations  is  Kriging.  Kriging  has  been  used  to  reduce  the  variance  of  regional  amplitude  ratios  for  earthquake 
populations  (Fan  et  al.,  1999,  Rodgers  et  al.,  1999).  Ordinary  Kriging  provides  a  minimum-mean-square-error 
prediction  (of,  in  our  case,  a  regional  amplitude  ratio  measured  at  a  particular  station)  at  a  given  spatial  location  that 
is  an  unbiased  linear  combination  of  previously  measured  data  at  N  separate  locations  (Cressie,  1993).  In  Ordinary 
Kriging  the  prediction  at  an  unobserved  point  that  is  located  far  from  any  data  points  converges  to  the  mean  of  the 
available  data  at  that  station.  For  a  station  that  has  only  data  in  a  subregion  with  anomalous  propagation  properties,  a 
prediction  in  another  subregion  far  removed  from  this  data  may  be  incorrect.  In  such  a  subregion  we  wish  for  the 
prediction  to  converge  to  a  value  that  is  the  average  over  all  stations  worldwide  of  the  same  region  type.  To  achieve 
this  limiting  behavior,  we  generalize  the  Ordinary  Kriging  approach  using  Bayesian  methods. 
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Figure  2.  Uncorrected  Pn/Lg(6-8  Hz)  and  Pn/Sn(6-8  Hz)  versus  distance  in  km  for  all  REB  events  and 
explosions  worldwide  in  tectonic  and  stable  regions.  The  solid  curves  are  best-fit,  three-parameter 
estimates  of  the  distance  dependence  for  earthquakes. 

Suppose  we  have  N  data  values  at  locations,  s ; .....  s  (V  for  a  particular  regional  amplitude  ratio,  say  Pn/Sn(6-8  Hz), 
recorded  by  a  single  station.  Let  x  =  (xi,...,xN)'  be  an  /V-dimensional  data  vector  where  the  x-t ,  i  =  l,...,N ,  are 
distance-corrected  values  of  log(P/S),  where  the  distance  corrections  are  relative  to  the  worldwide  wide  average 
distance  dependence,  estimated  by  using  P/S  values  for  data  available  from  all  stations.  We  normalize  the  data  so 
that  the  distance-corrected  worldwide  average  is  zero. 

We  assume  that  data  at  a  given  location,  s,-,  can  be  modeled  by  a  random  variable  with  a  normal  distribution,  with 

2  2 

unknown  mean,  /i(- ,  and  fixed  variance  (to  be  determined),  ,  where  o~  is  the  same  at  all  locations,  that  is, 

Xi  ~  .  (2) 

We  also  assume  that  the  means  of  the  normal  distributions  at  each  location,  ,  are  also  normally  distributed 

9 

random  variables  with  zero  mean,  equal  variances,  cc  ,  (to  be  determined),  and  correlation  between  any  two  flj  and 
/j  j  dependent  only  on  the  distance  separating  the  two  locations,  that  is, 

Pij  =COTT(pi,/J. j)  =  f(A(Sj,S  j))  ,  (3) 

where  A(Sj,Sj)  is  the  epicentral  distance  between  locations  s ,•  and  s j.  The  function / must  be  a  positive  definite 

function  (function  of  positive  type),  ensuring  that  the  covariance  matrix  of  the  random  vector  of  means  is  positive 
definite.  The  form  of  this  function  is  usually  taken  to  be  a  known  positive  definite  function  with  unknown 
parameters  estimated  using  the  available  data.  We  generally  use 

/(A)  =  exp(-A/a) ,  (4) 

where  a,  the  correlation  length,  is  estimated  using  the  data. 
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We  wish  to  determine  an  optimal  predictor  for  the  mean  of  an  unobserved  value  located  at  s(l  given  the  N  data 
values,  x,  located  at  .  We  wish  to  determine  the  distribution,  p(Po  I  x) ,  of  the  unknown  mean,  p() ,  at  the 

unobserved  location,  s(l ,  given  the  data,  x.  The  expected  value  of  this  distribution  (i.e.,  the  posterior  Bayes 
estimator)  could  be  used  as  the  estimator  of  p 0  at  that  location.  Under  our  assumptions  we  know  p(x(hx  I  p0,i) 

and  p(/j. Q,i),  where  i  =  (/i  | . p  N )' ,  so  we  can  use  Bayes  theorem  to  calculate  p(p o  I  x) .  Performing  the 

appropriate  integrals  we  get 

P(Pq  lx)°=  j[j p(x0,x\  pQ,i)dx0]  p(p0,l)di  .  (5) 

Since  integrals  of  normal  distributions  are  normal  distributions,  pip 0  I  x)  will  also  be  normally  distributed,  so  the 
proportionality  constant  need  not  be  determined  explicitly.  Substituting  the  appropriate  distributions  for 
p(jc'0,x  and  p(p Q,i)  as  described  above  and  performing  the  integral  over  Xq  gives 


pip Olx)oc  J 


exp 


exp 

—  -i'r’i 

2 

d  i , 


(6) 


where  y  and  i  are  (AT-l)-dimensional  vectors  given  by  y  =  (0,x')'  and  ^  =  (JllQ,i,),,  is  an  (AT-l)-by-(AT-l)- 
dimensional  matrix  given  by 


A1=J2 


0  0 
0 


(7) 


where  IN  is  the  identity  matrix  in  N  dimensions,  and  Z  is  an  t/V+ 1  )-by-(Af+ 1  )-dimensional  covariance  matrix  such 
that 


Sy=^cPy.  ;,  j  =  0,1,...,  A/ . 

The  integral  in  equation  (6)  can  be  shown  to  be  proportional  to 


pip o  I  x)  J  exp 


d  i , 


where 


and 


S-1  =Z“1+A1 


i=sai7. 


(8) 

(9) 

(10) 

(11) 


This  integral  is  easily  computed  by  using  the  theorem  that  if  a  random  vector,  i ,  is  distributed  as  ATX. S’) ,  then  the 
marginal  distribution  of  any  set  of  components  of  i  is  multivariate  normal  with  means,  variances,  and  covariances 
obtained  by  taking  the  corresponding  components  of  I ,  and  S  .  This  gives 


1  —1  —2 
~  — (S  00  )  (p  o-Po) 


PiPo  1  x)  exP| 

or,  equivalently,  p(j ,  given  x,  is  distributed  as  the  normal,  n(pq.ct2). 


where 


Po  -  (s  ^iy)o 


(12) 


(13) 


and 


a2  =  5oo=[(S“1  +  Tl1)-1]oo- 


(14) 
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A  prediction  surface,  which  can  also  serve  as  a  correction  surface,  is  obtained  by  computing  /)q(Sq)  for  all  values 

9 

of  Sq  within  regional  distances  of  the  station.  An  uncertainty  surface  is  given  by  <7  (Sq)  .  In  principle,  such  surfaces 
can  be  obtained  for  both  earthquakes  and  explosions  at  each  station,  however,  due  to  the  scarcity  of  explosion  data, 
only  earthquake  surfaces  are  practical. 

An  explosion  data  value,  say  Pn/Sn(6-8  Hz),  at  location  Sq  ,  can  be  compared  to  the  earthquake  surface  by 
computing  the  difference  x(Sq)- P.q(Sq)  .  This  should  be  a  large  positive  number  if  explosions  are  to  be 
distinguished  from  earthquakes.  Any  unknown  event  at  s()  can  also  be  compared  to  the  surface  in  the  same  way, 
however,  if  we  wish  to  compare  any  of  the  N  earthquakes  used  to  determine  the  surface  with  a  prediction  at  that 
location,  we  must  remove  that  earthquake  from  the  data  set  and  compute  a  new  surface  value  at  that  location  using 
the  N- 1  remaining  earthquakes  (the  leave-one -out  procedure).  There  are  three  parameters  that  must  be  estimated 

using  all  available  data,  <7C  ,  the  calibration  variance,  crr  ,  the  residual  variance,  and  a,  the  correlation  length. 

2  2 

Preliminary  work  suggests  that  a  is  typically  about  6  degrees,  and  that  a~  and  <7r  are  about  equal,  their  sum  being 
equal  to  the  total  variance  of  the  data. 

Figure  3  shows  a  typical  correction  surface  at  station  NIL  for  distance-corrected  log[Pn/Sn(6-8  Hz)],  with  the 
corresponding  uncertainty  surface  shown  in  Figure  4.  At  locations  far  from  available  data,  the  surface  approaches 
the  distance-corrected  worldwide  average.  The  locations  of  nuclear  explosions  and  earthquakes  recorded  by  NIL  are 
depicted  by  triangles  and  crosses,  respectively.  White  markers  indicate  values  that  are  less  than  the  surface,  while 
black  markers  indicate  that  the  value  is  greater  than  the  surface.  The  size  of  the  marker  is  proportional  to  the 
absolute  value  of  the  distance  from  the  surface.  Thus,  very  explosion-like  values  correspond  to  large  black  triangles. 
Figure  4  shows  that  the  uncertainty  is  lower  in  areas  containing  many  data  points,  while  it  is  higher  in  areas  around 
NIL  that  are  relatively  aseismic. 
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Figure  3.  Bayesian  Kriged  surface  for  log[Pn/Sn(6-8  Hz)]  at  NIL.  Triangles  show  locations  of  nuclear 
explosions  and  crosses  show  the  locations  of  earthquakes.  Black  (white)  markers  have  values  above 
(below)  the  surface.  Markers  size  is  proportional  to  the  absolute  distance  from  the  surface. 
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Figure  4.  Uncertainty  (variance)  surface  for  Iog[Pn/Sn(6-8  Hz)]  at  NIL.  Triangles  show  locations  of  nuclear 
explosions  and  crosses  show  the  locations  of  earthquakes.  The  markers  are  the  same  as  Figure  3. 

Event-Screening  Criterion  and  Score 


An  event-screening  criterion  can  be  developed  from  a  hypothesis  test,  using  the  corrected  data  and  the  associated 
uncertainty,  that  has  a  fixed  significance  level  with  respect  to  screening  out  an  explosion.  For  each  station, 
correction  surfaces  are  calculated  for  both  log[Pn/Sn(6-8  Hz)]  and  log[Pn/Lg(6-8  Hz)],  as  well  as  uncertainty 
surfaces.  To  construct  a  hypothesis  test,  we  define  a  scaled  variable,  X,  given  by 


y~  Hex 

sja~  +<7“Eq 


(15) 


where 


y  =  x(s0)  -  ju0(s0) .  (16) 

In  Equations  (15)  and  (16),  x(s0)  is  the  distance-corrected  value  of  log[Pn/Smax(6-8  Hz)]  for  a  test  event  located  at 
Sq  ,  where  Pn/Smax  is  Pn/Sn  if  Sn  is  larger  than  Lg  or  Pn/Lg  otherwise,  /I0(Sq)  is  the  correction  surface  value  at 

Sq  for  Pn/Smax,  /t|.x  is  the  mean  value  of  y  for  all  explosions  in  the  data  set,  0  ~  is  the  uncertainty  surface  value  at 
2 

Sq,  and  ex  's  the  residual  variance  for  explosions.  If  all  assumptions  are  true,  then  X  should  be  normally 
distributed  with  mean  zero  and  variance  one.  With  this  in  mind,  we  define  the  screening  criterion  to  be 

X  <  -za  ,  (17) 

where  za  is  the  (l-a)-percentile  of  the  standard  normal  distribution,  which  we  usually  take  as  2.576  for  a  0.005 

significance  level.  We  also  define  a  regional  score  such  that  an  event  is  screened  out  at  the  a  significance  level  if  the 
score  is  greater  than  zero. 
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Score 


(18) 


We  have  applied  this  screening  criterion  to  140  explosions  and  4173  REB  events  above  mb  3.5.  For  each  explosion, 
the  Kriged  correction  and  uncertainty  surface  values  are  computed  at  the  explosion  location  for  each  station  with 
regional  recordings,  using  all  earthquakes  available  at  that  station.  The  value  y,  given  by  equation  (2),  is  then 
computed,  from  which  the  mean,  /i|.x ,  can  be  estimated.  The  explosion  data  is  also  used  to  estimate  the  explosion 

residual  variance,  <7r  EX  >  which  we  estimate  to  be  (0.22)~.  The  other  estimates  used  are  ac  =  0.25 ,  <7r  Eq  =  0.25  , 

and  a  =  6.0  degrees.  For  each  earthquake  a  Kriged  and  uncertainty  surface  is  computed  at  the  location  of  the 
earthquake  by  using  the  remaining  earthquakes  recorded  at  that  station  (using  the  leave-one -out  procedure).  The 
scaled  value  X  for  each  event  is  determined  and  a  score  for  each  event  is  computed.  Those  events  with  positive 
scores  are  screened  out. 


For  our  data  set,  none  of  the  140  explosions  are  screened  out  at  the  0.005  significance  level,  while  about  73%  of  the 
4173  REB  earthquakes  are  screened  out.  Figure  5  shows  a  plot  of  the  scores  for  all  of  these  earthquakes  and 
explosions.  For  comparison.  Figure  6  shows  a  plot  of  distance-corrected  log[Pn/Smax(6-8  Hz)]  for  the  same  events, 
before  the  Kriged  corrections  were  applied.  As  can  be  seen  from  Figure  6,  if  the  explosion  with  the  smallest 
amplitude  ratio  had  been  used  as  the  screening  threshold,  only  about  50%  of  the  earthquakes  would  have  been 
screened  out,  showing  the  improvement  resulting  from  Kriging. 
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Figure  5.  Histograms  of  regional  screening  scores  for  140  explosions  (left)  and  4173  earthquakes  (right).  None 
of  the  explosions  are  screened  out  (all  have  scores  less  than  zero  using  a  significance  level  of  0.005), 
while  73%  of  the  REB  earthquakes  are  screened  out  (scores  greater  than  zero). 
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Figure  6.  Histograms  of  Iog[Pn/Smax(6-8  Hz)]  with  worldwide  distance  corrections  for  140  explosions  (right) 
and  4173  earthquakes  (left).  If  the  explosion  with  the  smallest  value  were  used  as  a  threshold  only 
about  50%  of  the  earthquakes  would  be  screened  out. 


CONCLUSIONS  AND  RECOMMENDATIONS 


High-frequency,  regional  seismic  P/S  amplitude  ratios  have  potential  utility  to  augment  the  seismic  event-screening 
procedures  based  on  event  depth  and  Ms:mb.  Screening  performance  is  generally  improved  by  correcting  regional 
amplitude  ratios  using  spatial  techniques,  which  help  account  for  path  variations.  In  this  paper  we  developed  a 
Bayesian  Kriging  algorithm  that  treats  subregional  variations  and  accounts  for  local  uncertainty  in  such  a  way  that  in 
regions  far  from  any  calibration  data  the  correction  surface  approaches  the  worldwide  average  and  the  uncertainty 
surface  approaches  maximum  uncertainty.  Using  the  screening  procedure  described  in  this  paper,  no  available 
underground  explosions  were  screened  out  at  the  Nevada,  Lop  Nor,  Semipalatinsk,  Indian,  or  Novaya  Zemlya  test 
sites,  while  about  73%  of  REB  events  above  mb  3.5  and  with  Pn-SNR  >  2.0  and  S-SNR  >1.2  can  be  screened  out,  a 
great  improvement  over  using  worldwide  distance  corrections  alone.  This  screening  performance  compares 
favorably  with  that  of  Ms:mb.  The  experimental  depth  and  Ms:mb  screening  criteria  at  the  pIDC  currently  screen 
out  less  than  50%  of  REB  events  above  mb  3.5  and  about  70%  of  REB  events  that  have  applicable  regional  data  are 
relatively  shallow  and  do  not  have  an  Ms  measurement. 

There  are  several  potential  ways  to  further  improve  the  regional  screening  performance.  The  hypothesis  test 
presented  in  this  paper  uses  only  one  screening  parameter,  namely  Pn/Smax(6-8  Hz).  Extending  the  test  to  include 
additional  amplitude  ratios,  e.g.,  Pn/Smax(4-6  Hz)  and/or  Pn/Smax(8-10  Hz),  generally  improves  screening 
performance.  Another  method  that  may  improve  performance  is  multistation  averaging,  where  amplitudes  ratios  for 
the  same  event  recorded  at  more  than  one  station  are  averaged  before  applying  the  hypothesis  test.  Currently,  about 
25%  of  all  REB  events  with  regional  amplitudes  have  recordings  at  two  stations  or  more,  a  percentage  that  will 
increase  as  more  IMS  seismic  stations  come  online.  Currently,  the  pIDC  measures  regional  phase  amplitudes  only 
on  vertical  components.  Although  P/S  ratios  measured  on  vertical  components  show  reasonable  power  for  separating 
earthquakes  from  explosions,  Kim  et  al.  (1997)  showed  that  P/S  spectral  ratios  of  rotated,  three-component  (3-C) 
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records  improve  the  power.  Since  most  IMS  seismic  stations  have  3-C  broadband  sensors,  the  horizontal 
components  should  also  be  considered  for  event  characterization.  Preliminary  work  is  this  area  has  shown  promise. 


Key  Words:  CTBT  monitoring;  regional  seismic  data;  event  characterization;  underground  explosions;  Bayesian 
Kriging 
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